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1. INTRODUCTION 


This report documents results accomplished at the NASA, National Space 
Technology Laboratories, Earth Resources Laboratory (NSTL/ERL) from the devel- 
opment and testing of an agricultural land use change monitoring capability as 
a part of the Land Cover Information Systems (LCIS) task of the AgRISTARS 
Domestic Crops and Land Cover (DCLC) project in cooperation with the U.S. 
Department of Agriculture, Statistical Reporting Service { SRS ) . 

Change detection analyses using remotely sensed digital data have been 
f plied to a variety of natural resource problems. They have been used in the 
monitoring of alteration in coastal ^ 9 » 26 , forestland^ » 15, 18^ rangeland, 
desert 4 * 20 , and interior wetland 8 environments as well as the measurement of 
land use dynamics in both urban 9 * 22 » 28 » 84 and natural settings 8 * 

This Information has been used to monitor for water quality changes in water- 
sheds'* 4 * l 6 and for increases in strip-mined lands 1 * ^ . 

The remote detection of locational changes in surface cover materials pre- 
supposes there are associated, measurable radiometric differences between 
successive dates corresponding to these changes on the ground. Geometric 
relationships preserved by imaging scanners make this possible. As a result, 
several methods have been devised to recognize and map these phenomena. Those 
tested have included: band ratioing 8 * 9 » 24 , band subtraction (image dif- 
ferencing)^* ^ 8 » 22 , pre-classification differencing (delta data classifica- 
tion^, post classification comparisons, 8, 10, 15, 23, 26, classification of 
multi date data^’ and measurement of spectral change vecto K s^’ ^ ^ • 



Other comparisons for change have been accomplished through con/ersion of the 
digital counts to absolute physical quantities that are then subtracted-^ and 
also by taking statistical measures of correlation, covariance, and/or 
percent-expl ai ned-variance by the first eigenvalue to compare between data 
sets^. One Investigator used the Kolmogorov-Smi rnoff (K-S) test to identify 
changes between dates 6 , while yet another has correlated land use changes with 
information contained in the third principal component of a transformed, multi- 
date data set 3 . 

Numerous obstacles prevent the straightforward execution )f these opera- 
tions, and subsequent problems make evaluation difficult. Basically, any of 
these methods require spectral data sets to be precisely co-registered so that 
the radiometric response of corresponding ground areas can be compared. Posi- 
tional inaccuracy adversely affects performance, although one method employing 
the K-S test reports that it is relatively independent of small misi^gistra- 
tion errors^. Other problems include the influence of time-dependent varia- 
tions of the extrinsic factors listed in Table 1. These factors variably 
combine to alter the radiometric fidelity of the recorded spectral response of 
a scene. This degrades technique performance by inducing the detection of 
untargeted factors. A few of these (e.g., clouds, cover material, and soil 
moisture changes) are locational by nature whereas other changes affect total 
coverage. Investigators have experimented with and applied various data modi- 
fications in attempts to negate or compensate for such factors. Generally a 
great deal must be assumed, and only the major influences are tr°ated, usually 
by "standardizing" or equalizing the effects on each data set rather than 
"correcting" or removing it as a factor since no one has proven that the 
effects can be entirely subtracted. 


t 


Table 1. Considerations for Temporally Dependent Sources of Change 
Reflectance Between Data Sets 


In 


Atmospheric Differences 

Clouds 

Haze 

Humidity 

Dust 

Seasonal Differences 


Solar Illumination Angle 
Phenol ogic St^.ge 

Surface Differences 

Soil Moisture 
Cover Materials 

Sensors/Systems Differences 

Orbital Altitude 
Platform Attitude 

Differential System Deterioration Rates 
Sensor Calibration 


Processing Differences 
Formatti ng 

Resampling Procedures 
Decompression Procedures 

Astrophysical Differences 

Solar Flux 

Magnetospheric Interference 
Various Axial Motion Components 
Ecliptic Variations 
Eccentricities in Orbit 
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Geometric corrections for skew, rotation, scaling, etc., are regularly 
applied to rectify both data sets In early stages of preprocessing before 
co-reg1strat1on5, 6, 9, 12, 17, 18, 19, 20, 26. other Investigators have 
tried to correct for aerosols (haze) 12 and other non-local 1 zed atmospheric 
effects20 as well as cloods3, 12. Atmospheric correction models usually 
require additional measurements acq.ired concurrently with the spectral data; 
otherwise data sets are standardized to one another by the same factors that 
are Inherently collected in the data^. If not for the unavailability of 
these concurrent atmospheric measurements and the expense and complexity of 
mathematically describing atmospheric Interactions, the use of these models 
would be more prevalent. More frequently used are standardizations for 
effects of differences In solar 111 umi nation angle and of differences In 
sensor calibration^* ^0, ^1 , 25. jh e e ff ec * s 0 f these factors can be 
eliminated, however, by selecting data collected on anniversary dates with the 
same sensor to avoid the major Influences of seasonal and system differences. 
But there are no guarantees that such factors as phenologlc conditions or soil 
moisture are as cyclical as sun angle effects and consequently have equal 
influences in each date. Also, data within these constraints of geographic 
and temporal coincidence may not be available, However, the data were avail- 
able for this study, and overall scene characteristics of both data sets were 
very similar. The author achieved better results with this approach when 
testing the first two techniques described in Section 3 than with data sets 
standardized by either the ERIM-developed coefficients 12 or the Landsat User's 
Handbook^. 21, 25 coefficients for both sun angle and sensor gain calibra- 
tion. 
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Another standardization based on the same theory as the regression 
modeling technique described In Section 3.4 was also tested. Attempts were 
made to model no change areas between dates in order to relate second-date 
spectral response as a function of the first date's response. This 
relationship was then extended to tlv entire second-date scene to describe It 
In the same terms as the first, which serves to equalize the extrinsic effects 
mentioned before. The transformation using the model coefficients worked well 
In the forested areas where the manual sample selection was adequate. A 
suitable sample of other types of unchanged areas of the scene was absent, and 
the model could not sufficiently describe these areas. This data set- 
dependent approach to standardization hypothetically enables scenes, regard- 
less of differences, to be related under equivalent terms. 

In addition to the problems of registration, extrinsic factor effects, and 
acceptable digital data selection, a control area Is essential to establish 
performance levels and to verify results when operating in an experimental 
mode. This inexorably creates the need for ground data coincidental in time 
and space with the necessarily retrospective study Interval to be used as a 
standard for comparison. Therefore definition of a control area dictates 
common areal coverage at two points In time from two different data sources. 
This restriction has prevented many investigators from having adequately 
substantive proof of performance from their results. Most often they do not 
satisfy all these criteria because the data are simply non-existent, or they 
are more interested in the success of the application. Usually aircraft 
photography has been the only reliable, alternative data source for large area 
surveys of this nature on successive dates, but even then coverage has been 
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spatially limited and Irregularly collected. In spite of this, aerophoto- 

13 16 

graphy has been used successfully to map land cover and land use change ’ 
and previously has been the only other feasible recourse until the advent of 
this technology. The control area, with Its attendant data requirements, Is 
necessary only for the experimental process. It Is not required for the 
routine application of these methods. 


Schemes for detecting surface changes take two basically different direc- 
tions In approaching their performance objectives. One, a technique Is used 
to examine or sample the entire pixel population as an undifferentiated set. 
As a result It may Identify anything from relative measures based solely upon 
the radiometric count difference of matched resolvable elements at two points 
in time 4 , or the technique may be able to locate and quantify specific types 
of change areas which indicate conditions at each time frame c * . Two, a 
technique may operate on a specified subset or stratum of cells where the 
changes known to be occurring are the subject of study such as In a particular 
ecosystem or habitat, with all other pixels being eliminated from analysis. 
This simplifies the procedure, and less confusion develops at the outcome 
because of the reduced number of data points to analyze. This allows the 
detection operations to address more subtle differences than could be recog- 
nized otherwise where these differences might be masked by greater spectral 
differences, which may not be of concern for a particular application, 
occurring in the general pixel population. Of course a means of different- 
iating the population is required prior to operating on the correct subset. 
The choice of schemes is dependent upon the objectives of the analyst for the 
specific purpose to be undertaken and potential for its best results. 
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The methods that were developed and tested here were Intended to express 
these changes In terms of the naturally vegetated landscape which underwent a 
change of conditions associated with agricultural production. Thus a survey 
of the whole population of data cells was taken as in the first option men- 
tioned. The objectives of technique execution were to detect and to locate 
changed pixel areas as well as to describe conditions at each date capable of 
characterizing the changes present with the least amount of ground truth. 
This does not mean that other methods or approaches will not provide these 
same informational components: detection, location, and identification. It 
was also desirable that they be flexible enough to provide maps as well as 
tabular accounts of specific change types that would affect agricultural 
productivity estimation. 
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2.1 Study Site Description 

The test site covers an area from 91.5° to 92° W. longitude and from 31° 
to 32° N. latitude on the fertile alluvial plains of the lower Mississippi 
River in east central Louisiana. This area is further characterized by 
minimal relief, poor drainage, and fertile, shallow, undeveloped, organic-rich 
soils complexly distributed by the fluvial processes at work in this region. 
Extensive mixed-bottomland hardwood forests of oak, gum, and cypress once 
dominated the landscape. Rapid clearing for agricultural production of crops 
and livestock has left less than a quarter of these forests standing 2 . Table 2 
shows the deforestation rates for parishes that are a part of this 


Table 2. Past and Predicted Acreages of Bottomla. •' ’ardwoods in Those Parishes Partly or Wholly Contained 
in Study Area. 
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study area. The magnitude and type of land cover change occurring within this 

area played a significant role In Its selection as an exploratory site for 

technique development and testing. 

2 . 2 Landsat MSS data 

Landsat MSS data sets collected October 10, 1974, and October 2, 1979, 

were obtained covering the frame deMned by path and row coordinates 25/38 of 

the Worldwide Reference System for Landsats 1, 2, and 3. Fall data sets were 
used as this Is the driest time of the year for this locale. It had been 
hypothesized from previous study 10 that forest and agricultural land cover 
conditions In this season would offer enhanced spectral separability by 
minimizing uk> problem of spectral overlap partly caused by the excessive 
surface wetness that prevails most of the year. 

The 1974 data set is in the pre-EDIPS X format for CCT's {57m by 79m 

resolution cell, geometrically uncorrected) while the 1979 data set is in the 

EDIPS P format of partially corrected, 57m X 57m resolution cells. A tech- 

27 

nique for overlaying Landsat data with Seasat data described by Wu was used 
to co-registcr and merge these differing data format types into a single 8 
channel, multidate source file. Registration was accomplished to within one 
pixel (57m RMS) of the base set. 

2.3 Aerophotographic Dat a 

Conventional, high altitude, color IR photography was available for 
retrospective ground coverage contiguously defining the test site. Aero- 
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photographic missions had previously acquired these data October 4, 1974, at 
1:120,000 scale and October 24, 1979, at 1:60,000. This closely coincided 
with the endpoints of the 5-year interval between Laidsat overpasses when the 
spectral digital data were acquired. (See Table 3.) 

The two sets of photography were analyzed for changes in land cover dis- 
tributions. Changed areas were delineated and rectified upon a common base 
using the eight USGS 15' series of topographic maps that comprise the study 
area. Digital land use change maps were produced using an X-Y digitizer to 
define the polygonal boundaries in the UTM coordinate system. Next this 
polygonal information was converted into a raster data file with each data 
cell representing either a change or no change area. This sequence is illus- 
trated for a 15' subset of the area in Paragraph 4.1. 

Table 3. Data Acquisition Dates for Data Types Used in 
This Investigation 



DATE ACQUIRED 

DATA SOURCE 

TIME 1 

TIME 2 

LANDSAT MSS 

October 10, 1974 

October 2, 1979 

FALSE COLOR IR 



AEROPH0T0GRAFHY 

October 4, 19/4 

October 24 , 1 979 


2.4 Map Information for USDA/SRS Areal Unit Analysis 


Prior to sampling and estimating the crops, livestock, farm labor, etc., 
for a geographic area, SRS personnel stratify land uses and agricultural land 
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use intensities through aerophotographlc interpretation. These land use 
strata are further subdivided Into what is the basic SRS geographic analysis 
area or frame unit. Each frame unit is a similarly homogeneous part of the 
stratum it represents. It corresponds to a specific polygonal area of the 
land surface bounded by permanent, recognizable map features. 

Strata maps (the area sampling frames) for Catahoula and Concordia Parish- 
es were obtained, and the frame unit boundaries were digitized. Superimposing 
these areal reference units with information produced by one of the change 
detection processes gives land use change statistics, as well as land use 
proportions at both dates, for each frame unit. Also, the geographic location 
and extent can be mapped on peripheral devices and used effectively to update 
the stratification. (See Figure 1.) here the same boundaries are used, and 
frame units only are redefined, which prevents reconstructing the entire frame 
by present means. 


3. CHANGE MONITORING TECHNIQUES 
3.1 Post-Classification Comparison 


This approach is one of the most widely used. It involves making inde- 
pendent land cover classifications for both points in time, usually through 
automated spectral pattern recognition techniques. These are reduced to 
common descriptive categories and then compared for areas of each category 
that have changed during the period covered. 
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FIGURE (Ic) SHOWS HOW IANOSAT MSS DERIVED LAND USES WITHIN USDA/SRS POLYGON*’. GROUND 
REFERENCE DATA (lb). WERE USED TO UPDATE THE LAND USE INTENSITIES OF THE AREA SAMPLING FRAME 
(la) FOR CONCORDIA PARISH. LA 


NASA/NSIl 

EARTH RESOURCES LABORATORY 


Flgura 1. Update of USDA/SRS Frame Units tor Concordia Parish 
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In this case an unsupervised clustering technique that passes a 
user-deMned window through the data to find spectrally homogeneous areas was 
used 10 . These are reduced Into statistically defined spectral groupings or 
signatures which provide the decision boundaries for mapping all data cells 
into classes based upon probability densities via a maximum likelihood 
algorithm. The 1974 data set produced 49 spectral classes, and the 1979 data 
set produced 54 such groups. Through Interpretive examination both were 
reduced Into the three major surface cove s that exist at this site: cropland, 
forestland, and water. With both time periods commonly represented by this 
classification scheme, they were numerically recoded in order to be digitally 
compared for detecting the desired changes. Figure 2 graphically depicts the 
assignment of each spectral signature to one of the three major surface covers 
and the position of the means of each signature on a plot of a visible and 
infrared band. 

3.2 Spectral Change Pattern Analysis 

This method uses the same pattern re-ogni tion technique in 3.1 as the 
primary data reduction meti d; however, rather than operating on individual 
dates, the co-registered, composite, multidate file or subsets of correspond- 
ing channels in each date are used as the source for statistical signature 
development. In this manner, as in the sense of mul tiseasonal classifications, 
the added information from another date defines classes whose spectra have 
changed in a distinct pattern in addition to those groups that have the same 
spectral response in both dates. In this way these composite signatures can 
be temporally sliced to indicate the spectral response at each date because 
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Figure 2. Plots of Spectral Group Means 
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the same spectral class represents both dates,, unlike the previous technique 
in which each date had signatures developed separately. Thus after mapping 
the data into these statistically defined categories, areas of change and no 
change can be identified by their signature migration along with conditions 
indicative of the surface cover at each point in time. 

In this investigation, the full eight channels of data, a six-channel 
subset (bands 4, 5, 7 of both dates) and a four-channel subset (bands 5 and 7 
for both dates), were tested and developed sets of 56, 52, and 58 training 
statistics, respectively. The results from the three data sets differ by only 
0.3%, with the four-channel data set giving the most accurate results. 
Whether this is because of the high correlation between MSS bands providing 
essentially most of the information in two channels or of the optimization of 
the classifies for four channels is not certain. Figure 3 graphically shows 
the assignment and position of each of the 58 spectral group means in two 
bands for both dates of the four-channel data set. Note here that the same 
class occurs for each date whereas before (3.1) the classes between dates have 
no relationship other than they represent corresponding types of surface 
material s. 


3.3 Radiance Vector Shift 


This method uses an algorithm that looks at two channels of information 
from each date Independently and then compares for differences between data 
sets in terms of standard deviation and an angle of relative, directional 
shift^. The algorithm works by finding all the pairwise occurrences from 
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Figure 3. Plots of Composite Signature Means 
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the two channels that have been selected from the first data set within the 
second, and It computes a distance from these corresponding pixel locations in 
standard deviations for that set of specific pixels. It then takes an angular 
measure relative to the first point, to the position of the second point ^hich 
is described by its values In the two channels. It was believed In the design 
of this particular algorithm this second channel of information would have 
descriptive value to the type of change, but as is, this directional component 
is not relatable through any common reference frame such as the origin of the 
two axes describing the Euclidean space that the values occupy. That Is, a 
number of possible land covers could have the same value for the bidirectional 
shift as well as equivalent measures of magnitude and be entirely different 
types of changes at both beginning and end. Another shortcoming to this 
algorithm's treatment of spectral change vector analysis is the averaging of 
the co-occurrence values in the comparison data set. Because of this, 
identity of the values resulting from that specific change is lost. 

A continuous range of change values is output where zero represents no 
change. A threshold is decided upon based upon ground data since there is no 
outstanding data feature to delineate the change/no change boundary. Usually 
this is gradational; the boundary may cover a range of five values or more. 
This condition leads to commission and/or omission errors wherever the thres- 
hold is set. In this case the ground truth was used to obtain the value 
producing the best results. 
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3.4 Regression Mod el 


This method Involves the development of a mathematical model through a 
stepwise regression procedure between each date that relates the second date 

^2) spectral response for 1 r Jividua1 ground cells to those for the first date 
Cl-|) for each corresponding channel of Information. The model values 

predicted for Tg > as described by its best fit with J\ t are subtracted from 
the actual T 2 data to produce a digital file of residual errors for all pixel 
locations. Areas of land cover change coincide with the more anomalous values 
derived from the predictive model. A critical value is determined for the 
residual error values, and pixels assigned a change/no change status 
accordingly. In all trials between corresponding bands in each date, the 

relationship was best described by a cubic equation in the form: 

v 2 

Y i j k = A o + A l^ijk + A 2 X ijk + A 3X J ijk + E 1jk 

where; 

Yjjk = band k value at row i and column j in second date (T 2 ), 

X-jjk = band k value at row i and column j in first date (T]K 
A 0 = constant offset, 

a 1 ,A 2 ,A 3 = multiplicative factor for first, second and third order 
regression coefficients, and 

= Y observed - Y predicted, residual error tha.c represents 
change to some degree beyond predicted fit between dates for 
the . ound area imaged at row i and column j in band k. 

The basic precept here is that if there were no change, E-jj^ = This would 
be the case if an area could be imaged twice in short succession before any 
measureable changes could take place or if this relationship was established 




for a duplicated data set: a (or T 2 ° T-j). But as the Interval 
between successively collected data Increases, this relationship evolves to 
express whatever changed conditions present can be mathematically described -- 
In this case, a cubic expression. 


Earlier trials did not deliver the expected results. Sampling of corres- 
ponding cells between dates had depended upo'i a coarse, regular Interval of 
point selection, because of program llmltiatlons and study area size, and 
proved to be Inadequate to describe the desired relationship. After reconsid- 
eration, it was decided that to properly describe T 2 response as a function of 
, the model should express the relationship between dates In terms of no 
change. In other words, the samples use.. to develop the model should be 
selected from areas with absolutely no location-specific surface changes. In 
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could be taken into consideration and be expressed by the model as a constant 
offset, and (2) the calculated, predicted values of the model would reflect no 
surface changes so that (3) In computing the residual errors between the 
model's predicted value and those actually observed, highly anomalous values 
would occur In areas of change. Unfortunately, only forested areas of no 

change could be stratified for model building and did not represent the 
entirety of land covers within the scene The manual or supervised procedure 
of sample selection for defining the model was unaMe to locate enough 
acceptably unchanged examples of other representative land covers to success- 
fully describe them through the stepwise regression analysis. Of the major 
land covers characterizing this area, no agricultural samples could be 
used--even though there were large areas of agriculture in both dates--because 
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of the continual changes in surface conditions brought about by their 

intensive human use. Even the large areas of water in this area were 
rejected, because they also varied extremely between dates. With no other 
type of samples but forest to develop the relation between dates, the model 

did not perform well under these other condtions. However, within forested 

areas, the method worked tsonably well. There was very little error or 

noise, and with experience the various residual error levels could be 
associated with specific types of change. 

In expressing the locational, no-change relationship between dates, this 
method might also work well as a data-specific standardization between any 
data sets from which the model was developed. TMs idea was not fully tested 
because of the same sampling problem. 

4. EVALUATION OF RESULTS 
4.1 Verification Procedure 


For verification of detected changes of land use between 1974 and 1979, 
the digital ground data mapped from the coincident aerophoto coverage was 
formatted into a multichannel, georeferenced data file. Every 57m^ cell 
within polygon boundaries representing the photointerpreted land use change 
was encoded with the value "1". All areas of no land use change were assigned 
a zero value. This exercise served to provide complete, contiguous data 
representation for an area of eight 15' series quadrangle maps with a digital 
land use change map to serve as a comparison standard for method performance. 
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Each Individual change detection technique's digital output was also reg- 
istered into the data base as "0's" or "l's". This was accomplished after 

7 

registering the Landsat data to the UTM coordinate system so that equivalent 
points on the ground could be compared. 

Criteria for photointerpretation employed the use of a minimum ten percent 
crown closure to distinguish forests from non-forests. Only change areas 
larger than ten acres were delineated during photointerpretation. In order to 
maintain as much label definition consistency and compatibility between data 
sources as possible, Landsat-derived data products were further subjected to a 
spatial classifier recognizing only change parcels larger than ten contiguous 
acres, effectively eliminating anything less from the comparison. Each tech- 
nique's output was then added cell-by-cell to that of the doubled ground data 
value to produce an "error source map" and/or accuracy statistics. This 
operation is better described by CHO ~ CHl + (CH2 * 2 ] 
where: CHO = result of operation for comparison 
CHI a MSS-derived change data 
CH2 = ground reference date 

for each equivalent ground resolution cell. Trr- possible outcome of this 
operation is a 0, 1, 2, or 3 for each cell where (0) zero represents agreement 
between both data sources that no change has occurred for that cell, (1) one 
indicates commission error on the part of the computer-identified change 
technique, (2) two indicates omission error, and (3) three indicates there was 
mutual agreement to that data cell having undergone changes. Statistical 
information from this procedure for the four techniques tested is shown in 
Table 4, and an "error source" map is shown in Figure 4. 


Table 4. Accuracy Tabulation of Verification Procedure ’Results 
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Figure 4. Accuracy Verification Procedure Between LandeatDerlvvd 
end Aerophoto-Derlved Lend Cover Change 
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Usfng the error-source map generated from this procedure, the correspond- 
ing locations of the commission and omission errors can be examined in the 
photography and a descriptive determination of their cause can be identified, 
whereas the numerical procedure alone only classifies the errors in one of 
these two ways. With this ability, the first examination revealed that many 
of the larger areas of commission error were locations of actual change that 
had been overlooked during the photointerpretation. The data in these were 
locations updated to correct these errors in the ground truth. After further 
reexamination of this information in its spatial context, the remaining errors 
were attributed to the following factors: 

1. Non-simul taneous acquisition of Landsat and aerophotographic data as 

in the case of sizable random, locational errors of omission; 

2. Several types of misregistration of the ground data to a common map 

base such as: 

(a) photo-to-map transfer of land use change delineations, 

(b) imprecise digitizing of these locations, 

(c) conversion of this polygonal data to a raster data file, as in 
small, contiguous errors of both commission and omission in 
boundary locations; 

3. Misregistration between data types such as: 

(a) band-to-band registration in individual spectral data sets, 

(b) scene-to-scene registration between spectral data sets, 

(c) scene-to-map registration for a georeferenced data set which 
resulted in more scattered, but patterned, errors of both types in 
many physical boundary locations; 




4. Human error In ground data set development such as: 

(a) incomplete Identification of all change sites between sets of 
photography that could not be corrected and 

(b) misinterpretation of land use change which resulted in small 
discrete errors of both types; 

5. Discordant labeling criteria between data types--This factor caused 
either commission or omission errors depending on the circumstance. 
For instance, a computer-assisted, satellite-detected spectral change 
may consistently occur at a forest density break of 40 percent crown 
closure, whereas manual mapping criteria may stipulate a 10, 15, or 20 
percent break before it is recognized as a change to another category. 
Other problems of this nature included surface areas covered by high 
water and the range of surface conditions associated with cultivated 
areas. 

6. Spectral similarity between certain surface materials and consequent 
co-classification — Small examples of various misclassifications were 
found that included confusion between burned over areas and wet areas, 
between some types of agriculture and forested sloughs, brakes, and 
wetlands, and within highly complex boundary areas where many land 
cover types occur within a single resolution cell and produce 
integrated spectra. 

Most of these errors (1, 2, 4, and 5) could be eliminated in routine applica- 
tions where the verification exercise is unnecessary. 

4.2 Summary Conclusions 

Computer generated classifications of Landsat mul tispectral (MSS) data can 
be used to measure forestland to agricultural land use changes accurately 
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when the proper data are selected and land parcels being converted are ten 
acres or larger. With these stipulations, the accuracy obtainable Is at least 
equivalent to what can be obtained through aeropi.otographic measures of 
changed land use. Results showed that approximately 10.5% to 13% ( , the 
entire half degree by one degree study area had changed from forestland to 
agricultural use In the five-year period. However, in some areas where land 
use changes were more concentrated, data from certain 15' quadrangles 
Indicated up to 20% of this land use change within their boundaries. Several 
blocks of land as large as four square miles incurred 100% clearing and 
replacement by agricultural use. 

The methods reported here work without extensive efforts to standardize 
various extrinsic effects on each data set. All techniques require accurate 
digital co-registration of the data sets. The two methods involving maximum 
likelihood classifications as the primary dita reduction tool provide all the 
information requirements discussed in Section 1 with the least ground truth. 
These methods also more accurately depict the geograohic distribution of the 
changes identified. This is paramount in applications where this spatial 
detail is necessary. Even though all techniques' performance results vary 
from 89% to 96% correct and appear adequate, there is a marked difference in 
the images each technique's accuracy produces. (See Figure 5.) The added 
dimension of this additional information suggests that many applications would 
be unusable without at least a 95% accuracy by this method of accuracy 
measurement. 

The post-classification comparison (PCC) and spectral change pattern 
analysis (SCPA) techniques obtained the same scores for accurately identified 
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change and equal amounts of omission error because of similarities in the data 
reduction techniques used in each; however, the SCPA technique was more 
sensitive to decreases in forest cover density despite scoring less commission 
errors and despite user-supplied labels on spectral groups in the PCC tech- 
nique that eliminated differing labeling criteria. Reasons for this will be 
explored as these methods are tested in test areas in Kansas and Arizona. 
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